GAMMA-EXPANSION FOR A ID CONFINED LENNARD-JONES 
MODEL WITH POINT DEFECT 

THOMAS HUDSON 
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Abstract. We compute a rigorous asymptotic expansion of the energy of a point defect 
in a ID chain of atoms with second neighbour interactions. We propose the Confined 
Lennard- Jones model for interatomic interactions, where it is assumed that nearest neigh- 
bour potentials are globally convex and second neighbour potentials are globally concave. 
We derive the T-limit for the energy functional as the number of atoms per period tends 
to infinity and derive an explicit form for the first order term in a T-expansion in terms 
of an infinite cell problem. We prove exponential decay properties for minimisers of the 
energy in the infinite cell problem, suggesting that the perturbation to the deformation 
introduced by the defect is confined to a thin boundary layer. 



1. Introduction 

The analysis of discrete lattice systems and their relationship to continuum mechanics 
is currently a growing area of study within applied analysis. Many rigorous results have 
been obtained in the past ten years connecting discrete models with continuum limits, 
which are extensively surveyed in 



BLBL07 . The most well-developed approaches have 



been either to apply r-convergencdj to discrete energy functionals parametrised by the 
number of atoms per unit volume in the model (see Sch06 AC04 , BG02] ) , or to apply 
forms of the inverse function theorem to show that for a discrete energy with the same 
parameter fixed, the Cauchy-Born rule holds; i.e. for a given atomistic deformation and 
a certain range of atomic densities there exist continuum deformations which are close in 
some norm, and have a similar energy (see OT12,EM07|). 

Here, we take the former approach. We build upon recent works on surface energies in 



discrete systems SSZ11,BC07 , which employ T-development' as first defined in AB93 



and extensively discussed in BT08 . We define energy functionals with and without 
defects, and present the Confined Lennard- Jones model for interatomic interactions, which 
we motivate with a formal analysis. We then investigate the scaling of the perturbation to 
the energy which is introduced by the defect. We also provide a concrete cell problem that 
may be used for explicit computation of the first-order energy, and show that a minimiser 
of this cell problem decays exponentially away from the defect. This allows us to conclude 
that minimisers of the energies with and without a defect are essentially the same except 
on a thin boundary layer around the defect. 

1.1. Motivation. The tools developed to study the relationship between atomistic and 
continuum models rely upon the high level of symmetry which is maintained after deform- 
ing a crystal. However, the pure lattice behaviour is not the only factor in determining 
the bulk properties of such materials. The last century saw a revolution in the materials 
science community, as it was realised that lattice defects can change the strength of an 
otherwise perfect crystal by orders of magnitude. Understanding defects, how they scale 
and in what rigorous ways one might modify the continuum approximation of crystalline 
solids to take them into account is therefore key to developing our understanding of how 
best to model and predict their behaviour. 



For an introduction to T-convergence, see 



Bra02 
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(a) 2D Defects (b) ID Defect 

Figure 1. Impurities and Interstitials. 



As a first step towards this goal, we consider arguably the simplest crystalline defect, 
a dilute point defect. A point defect is an interruption of the pure lattice structure 
caused by changing an atom at some lattice site (called an impurity) or by inserting an 
atom into the structure at a point which is not a lattice site (called an interstitial). In 
ID, impurities and interstitials are essentially identical when atoms are treated as point 
particles with hard-core interactions, since atoms cannot move past one another, and so it 
is easy to modify the reference configuration to take such a defect into account. In higher 
dimensions, the two defects are qualitatively different (see Figure II]), with an interstitial 
requiring an additional point in the reference configuration which breaks the symmetry. 

The model we analyse is one-dimensional, and to avoid surface effects, we use a periodic 
reference domain, so that in effect the atoms lie on a one-dimensional torus. The defect 
considered is dilute since only one atom in the chain is of a different type. 

By computing the r-limit of the sequence of energy functionals for the model described 
here, we arrive at an energy which encodes some of the properties of the minimisers of the 
functionals along the sequence, but no quantitative information about the error made. 
Computing higher-order limits gives further, more quantitative control on the energy 
minima, and in this case will also allow us to say something more qualitative about the 
minimisers. 

1.2. Outline. As discussed above, we apply T-convergence to a sequence of atomistic 
energy functionals that depend on the parameter e, which is the inverse of the number of 
atoms per unit volume. 

In the remainder of Section [TJ we propose a model for interatomic interactions in a ID 
chain including a point defect. We then make a formal analysis of the model to motivate 
this study and the Confined Lennard- Jones model which we propose, before reformulating 
the problem in a format amenable to analysis in the framework of the one-dimensional 
Calculus of Variations. 

In Section [2j we derive the r-limit for the series of functionals defined in Section [I] as 
e tends to zero, and note that the introduction of the defect does not perturb the r-limit 
at this order. 

In Section [3j we collect and prove some results about the minimum problem for the 
O th -order r-limit of the atomistic energies, including existence, uniqueness and regularity 
for minimisers. 

In Section |4j we proceed to derive a first-order r-limit, expressing it in terms of a 
minimisation problem in an infinite cell, and prove some properties of this minimum 
problem along the way. 
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1.3. Physical model. For now, fix N G N. We consider 2N atoms indexed by 

ie{-N,...,N-l} 

which have a spatial density of 2N/L, where L is the total length of the deformed config- 
uration. To make the domain periodic, we identify an atom indexed by i = N with the 
atom i = —N so that we avoid boundary effects, and defining e = 1/2 AT, we choose to 
take reference positions for these atoms to be 

X{ . IE c Si e . ^ ^! 2 ' ' ' ' ' 2 J ' 

We also define 

n:=[-i|]» 

so that Q e = Qr\eZj. It should be noted at the outset that the choice to use 2 A" atoms will 
not be restrictive to our analysis but will make some of the concepts easier to elucidate, 
and that we will frequently write s — y to mean iV — y oo. 

Fixing the coordinate system so that atom —N lies at 0, any configuration can be 
described by a map 

y : tt £ -> [0, L] such that y{-\) = 0, y{\) = L. 

We will use the shorthand 

Vi ■= y(xi), 

and extend y to a map on the whole of eL by defining 

V2kN+i+l — V2kN+i '■= Vi+1 ~ Vi 

for any fceZ. 

The atoms in our model are assumed to interact through pair potentials which decay 
rapidly so that it suffices to consider an interaction between atoms and their 2 immediate 



neighbours on either side. As explained in Section 1.1 all atoms except one are of the 



same type, regarded as the 'pure' species. As in BDMG99 , the potential energy of a 

led tc 

Djyi 



bond between atoms is assumed to be expressed as a function of the relative displacement 

Ui+j ~ Vi _ Ui+j ~ Vi 



In the case where all atoms are of the pure species, a bond with relative length s has 
energy (j)i(t) for nearest neighbours, and fait) for second neighbours. The internal energy 
of the configuration arising from the interatomic forces is 

JV-l 



E £ p iy):= Yl MDiVi) + <h{I>2Vi 



-N 



Since in each configuration we assume there is a single defect, we assume without loss 
of generality that the defect is at index i = 0. The energy of bonds of relative length s 
between this atom and its neighbours are ij)\(t) for the nearest neighbours and ip2it) for 
second neighbours (see Figure [2]). The introduction of the defect causes a modification of 
the energy which is given by the addition of the following energy term: 

E d (y) := faiD 2 y- 2 ) - faiD 2 y- 2 ) + ^liDiV-x) - 0i(-Di3/-i) 
+V>i(-Dij/o) ~ 0i(-°iyo) + faiD 2 yo) - faiD 2 y ). 
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Figure 2. Pair Potentials 

Finally, we also consider dead loads fi acting on each atom. Taking as initial positions 
the points L(xi + |), the work done by these forces is 



JV-l 



E!(y):= Y, fi{Vi-L(*i + \))- 



This term can be though of as the work done by a linearisation of some external force 
field near the homogeneous linear state yi = L(x, + |). To keep notation concise, we will 
frequently write u, to mean 

Ui := yi -L(xi+ §), 

and we extend / by periodicity to a map over eL by defining 

J2kN+i — fi 

for any fcgZ. 

The total energy for the atomistic system considered is therefore 

E%y):=E;(y)+E d {y)+Ef(y). 

1.4. Formal analysis. We expect that atoms should minimise the energy E e , and we 
therefore seek to characterise the minimal energy and the states which attain this min- 
imum. We will consider the situation when N is large, and when the material is behaving 
elastically. In this case, interatomic displacements should vary slowly over the domain, 
and so we assume the Cauchy-Born hypothesis holds (for more information, see Zan96|). 



This states that interatomic displacements follow linear deformations of small volumes of 
the solid, and so we assume that 

D 1 y i ,D 2 y i ~ Dyfa), 

where y : Q — > [0,L] is some suitably smooth function describing the displacement. This 
means that the energy 

M D iVi) + M D 2Vi)-<i>i( D y( x i)) + M D y( x i))- 

Motivated by this we define W, the continuum elastic energy density to be 

W(t):=Mt) + Mt)- (1-1) 

The total energy is now approximately 

N-l 

E% y ) ~ Yl w ( D y&)) + E M + Ef(y)- 

i=-N 

We expect this energy to grow linearly as the number of atoms increases, so it makes 
sense to look at the mean energy per atom, eE £ , as N gets large. The size of the defect 
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is fixed and small, so eE d (y) should vanish as e — > 0, and 



eE%y)~e ]T W(Dy(x i ))+f i 



Mv)'= [ W(Dy) + fudx, 



where u(x) = y(x) — L(x + \). The minimiser of the right hand side should satisfy the 
Euler-Lagrange equation for this functional, 

A.(W'(Dy))=f(x). 

This equation can be integrated to give 

W'(Dy) = a(x) + E, where a(x) = [ f{t) dt. 

.7-1/2 

The defect should then contribute a further term proportional to its size, 0(e), to the 
energy. Since the mean energy is only perturbed by a small amount, we should expect 
that any perturbation to the minimiser would also occur close to the defect, due to the 
mismatch between 0j and ipi there. The defect always remains at % — 0, so when N is 
large, we make the ansatz that close to the defect D\yi ~ F Q + r i) where F := Dy(0) and 
Ti is a small perturbation. Defining 

o-ljv := -\ef-N and of := a\_ x + &U-\, 
then integrating by parts, we can rewrite the external force terms as 

N-l N-l 

y fiUi= y. a l D\Ui and / fudx = / aDudx. 

i=-N i=-N J Q ^ n 

This then means the additional contribution is 

£E£{y) ~ My) ^E 00 (r):=E d (r)+ J2 (&(F + r< ) + 2 (F O + ^±l) - W(F ) 

l = — OD 

+ Oi (F + n-L)- o-(O) (F - L)) , 
where we have defined 

E d (r) := MFo + ! -^ i ) - M F o + !I -^ i ) + M F o + r_i) - 0i(F o + r_i) 

+ 0!(F O + r ) - Vi(F + r ) + ^ 2 (F + ra±n) - 2 (F O + a±n). 

If / is smooth enough, then for % close to and e small, a\ ~ cr(0), so the integrated 
Euler-Lagrange equation for T§ gives 

oo 

£ 00 (r)~£ d (r) + J2 M F o + ^)+M F o+ r - 1± Y ±1 )- W ( F o)-{W'(F )-^)r t . 

i=—oo 

The sum of the E r* terms should vanish due to the boundary conditions. Since we expect 
decaying solutions, we linearise in r^ away from the defect, giving 

^(r) ~ F c (y) + \ J2 <fi(F )rj + 0' 2 '(Fo)(^) 2 . 

\i\>R 

For a minimiser of this energy, the r* should approximately satisfy 

i^(Fo)(ri_i + r<) + ^'(Fo)ri + |0 2 '(F o )(r, + r i+1 ) = 0. 
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Making the usual ansatz n = a A 4 , a solution must satisfy 

§0' 2 '(F O ) + W"(F )\ + i0' 2 '(F o )A 2 = 0. 

If 0i and W are convex and 02 is concave at Fq, as is the case in Lennard- Jones type pair 
potentials, then a straightforward analysis of the roots of this equation implies that there 
are two positive real roots which multiply to give 1. These two roots correspond to an 
exponentially decaying solution and an exponentially growing solution. 

Rigorous versions of these formal results will be the subject of this paper, and motivate 
the assumptions we make about the potentials and external force in the following section. 

1.5. Confined Lennard-Jones Model. Motivated by the formal analysis carried out 
in the previous section, this section details the assumptions that we make about the pair 
potentials and external force field. 

We will assume that all potentials 0j and ipi are C 2 on the interval (0, <x>). Additionally, 
we assume the potentials and external forces satisfy the following conditions. 

(1) The nearest neighbour potentials are infinite for negative bond lengths and blow 
up as bond lengths approach zero, i.e. 

0i(t),^i(t) = +oo for s < 0, 
lim0i(t) = +oo and lini'0i(t) = +oo. 

(2) The nearest neighbour potentials are /-convex, i.e. for any s > 0, 

, 1 '(t),<(t)>Z>O. 

(3) The second neighbour potentials 02 and ip2 are concave. 

(4) The second neighbour potentials 02 and -02 are 'dominated' by the nearest neigh- 
bour potentials 0i and ipi, i.e. there exist constants a £ (0, 1) and C £ M such 
that 

4>2{x) > -«0i(x) + C, <hi x ) > —onpi(x) + C, 

02(^) > — ol(J>i(x) + C, ip2(x) > -aipi(x) + C. 

(5) The 'pure' potentials are such that the resulting continuum elastic potential is 
/-convex, i.e. defining W as in (1.1[), for any t > 0, 



W"{t) = <p'[{t) + ^{t)>l. 
(6) We assume fa = f(xi) where / £ C 2 (Q). 

Remark 1. Assumption (1) prevents atoms from exchanging positions with respect to 
the reference configuration, and ensures that we prevent plastic deformation. 

Assumptions (2) and (3) are made to simulate the behaviour of a Lennard-Jones type 
potential, which is convex for short bond lengths, and then concave after for any bond 
length past some critical length; see for example Figure [3] When under strains in the 
elastic regime, the nearest neighbours lie in the convex part of the potential, and all other 
atoms lie in the concave part. 

Assumption (4) enforces the elastic behaviour of the material and prevents fracture 
from being favourable. 

Assumption (5) prevents any form of microstructure from forming, but since the decay 
of Lennard-Jones potentials used in applications is always relatively rapid, this assumption 
is reasonable, and for the sake of clarity we avoid significant complications to our analysis. 

□ 
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Figure 3. Possible choices of 0i and <p2 with the assumptions prescribed 
which approximate a Lennard- Jones potential, 0lj- 

These assumptions lead to the following facts which we will use frequently throughout 
this paper. The fact that <p2 is concave implies that 



iMa) + 1Mb) + 2 (^) > \W{a) + \W(b). 
By using the concavity of the second neighbour potentials again, we have 



:i.2) 



\M°) + M a -f) + Mb) + M b -^) + Mc) + M 9 ?) + \Md) 

> l{Ma)+Ma)) + {Mb) + 1Mb) + ^2(6)) 

+ {Mc) + \Mc) + |0 2 (c)) + \{Md) + Md)). 

Assumption (4), that the behaviour of nearest neighbour potentials dominates, implies 
that 

\M°) + M a -f) + Mb) + M b -¥) + Mc) + M c -f) + \Md) 

> (1 - a) (iMa) + Mb) + Mc) + \Md)) + 3C, 

> \l (|a 2 + b 2 + c 2 + \d 2 ^j + 3C, (1.3) 

where on the last line we have adjusted the definition of I to keep estimates concise 
throughout this paper. This estimate will allow us to prove coercivity results which 
ensure that sequences of deformations with uniformly bounded energies are compact. 

1.6. Function spaces and topologies. In this section we define the topologies with 
respect to which we will carry out our analysis. 

Throughout this paper, we use the usual notation for Lebesgue and Sobolev spaces, and 
use || • || p and || ■ || lp to denote the usual norms on L P (Q) and W 1,p (f2). We write y £ — )■ y 
in L p or y e — > y in W 1 ^ (or H 1 ) to mean 

\\y £ -y\\ p ^Q or \\y £ -y\\ hp -»■ 
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respectively. We will also refer to convergence in the weak topology on H 1 (fi); we say 
Ve ~~ ^ V or Ve converges weakly to y in H 1 if for any / e H _1 (f2), 

(f,y £ )^(f,y), 

where ( • , • ) is the usual inner product on H 1 (fi). 

Since we wish to find the r-limit of the sequence of energy functionals eE e defined 
above, we need to ensure they are defined over the same space. Following Braides, Dal 



Maso and Garroni in BDMG99 , we associate any discrete deformation y e with a piecewise 



linear interpolant defined everywhere on Q,, 

y £ (x) := y\ + D x y\ (x - x t ) for any x € {x h x i+1 ). 

For each choice of e, these linear interpolants lie in the spaces 

Pg(fi) := {y e W 1 '°°(Q) : y is linear on (xi,x i+1 ), Xi,x i+ i G Q £ }- 

The admissible deformations A e (L) are defined to be the set of such interpolants with the 
correct boundary conditions: 

A%L):={yePl(n):y(-l) = 0,y^) = L}. 

For any e, 

A £ (L) C A(L) := {y G H 1 ^) : y{-\) = 0,y(|) = L}; 

in fact, in the weak topology on H 1 (fi), it is well-known that the sequential closure 



\JML)=A{L). 



N=l 



In Section 2.1 , we will show that this topology arises from the assumptions made in 
Section 11.51 

We are now in a position to suitably extend the functionals eE e so that we can take 
a r-limit. Guided by other work for similar models (amongst others, see those used 
in 



BC07l|BG02[|BDMG9 9|) we define 7 e : A(L) -> R to be 

I +oo otherwise. 



2. The O th -order T-limit 
Our first result gives the first term in the T-expansion of T s 



Theorem 1 (O th -order r-limit). With respect to convergence in L 2 , 

r-lim F{y) = My) ■= f W(Dy) +f( y -L(x+ §)) dx. 

e ^ u Jn 

The T-convergence of the internal energy in this result is already covered by Theorem 3.1 



in BC07 , but we present a complete proof here in order to demonstrate the special 
structure of the Confined Lennard- Jones model described in Section L5 By exploiting 
the convexity and concavity of the potentials, we do not have to resort to a homogenisation 
formula to prove the liminf inequality. 

As with any T-convergence result, we need to prove the relevant liminf and limsup 
inequalities. We present the proofs of these inequalities in turn. 
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2.1. The liminf inequality. The liminf inequality is the following statement: 

Proposition 2. If y £ —¥ y in L 2 then 

hminfj-% e )> T (y). 

To prove this, we use the fact that if 

supj-' £ (y £ ) < +oo and y £ — > y in L 2 , 

then y £ — ^ y. This equicoercivity result is encoded in Lemma [3] We then prove that 
J-" £ (l/ £ ) is approximately bounded below by J^oiy 6 ) for any given deformation y £ e A £ (L), 
and finally we can use the fact that Tq is lower semicontinuous with respect to weak 
convergence in H 1 to obtain the inequality required. 

Lemma 3 (Weak H 1 coercivity). If y £ — > y in L 2 and ^(y 6 ) is uniformly bounded 
for all e > 0, then y £ — *■ y in H 1 . 

The proof of this result relies upon the growth assumptions and estimates made in 



Section 1.5, and is inspired by the argument used in the proof of Theorem 4.5 in Bra02 
Proof. First, we estimate the energy below away from the defect. For ease of reading, let 

P e :={-N,...,N-1}\ {-2,-1,0}, 
i.e. the set of indices which have only pure interactions with their two neighbours on the 



right. We now use the inequalities from the end of Section 1.5 The estimate made in 
(|1.2[), and the /-convexity of W imply that for any y E A £ (L) 



e Y, iMDiVi) + MDiVi) + -MDiVi+i) > e ^ \W{D lVi ) + \W{D iyi+x ), 

ieP £ 

>eJ2¥ (\ D M 2 + \Diy i+ i\ 2 ) + C. (2.1) 



i&Pe l&Pe 



ieP e 



The estimate made in (1.3) then allows us to bound the energy coming from bonds near 
the defect below. 



lMDw-2) + M^y-2) + MDw-i) + <M£>2y-i) + M D m) 



+ M^Vo) + \<k(PiVi) >J2l l (l^l 2 + l^i^+il 2 ) + C - ( 2 - 2 ) 

i=-2 

Combining these estimates, we have that 

JV-l 

eE £ p (y £ ) + eE d (y £ ) > e J^ (* \ D ^t\ 2 + d] 



i=-N 



l\\Dy £ \\l + C. (2.3) 



Finally, Lemma 5.3 in BDMG99| implies that 



eEf{y £ )^ [ f( y -L(x+±))dx, 



so it follows that eEf(y £ ) is uniformly bounded. Combining this fact with estimate (2.3) 
the uniform bound on T £ implies 

sup ||-D?/ e || 2 < +oo. 

£>0 
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The argument is now concluded via the standard result that y £ — *■ y in H 1 if and only if 
H^lli 2 is uniformly bounded and y £ — > y in L 2 (f2). D 

Remark 2. Lemma [3] can be interpreted as saying that if the mean energy is bounded 
along some sequence of atomistic deformations, then the interatomic strains do not get 
too large, since they are compact in the weak topology on H 1 (fi). This reinforces the 
notion that we are in an elastic regime. □ 

Lemma [3] permits us to use the fact that as W is /-convex, the map 

y\ — ► / W(Dy) dx (2.4) 



jci 

is lower semicontinuous with respect to weak convergence in H 1 (fi) (see for example 
Corollary 2.31 in Bra02|). Fix 5 > 0. The estimates made in (2.1) and (2.2) imply that 



sE £ p (y £ ) + sE d (y £ ) > / W(Dy £ ) dx + C 5. 

Ja\(-s,s) 



Using Lemma [3] and the weak lower semicontinuity of (|2.4), we have that 



lim inf (eE £ p (y £ ) + eE d {y £ )) > lim inf f W(Dy £ ) dx + C5 

£ ^° £ ^° JQ\(-S,S) 

> f W(Dy)dx + C5. 

Since 5 was arbitrary, we let 5 — > 0, giving 



lim inf (eE £ p (y £ ) + eE d {y £ )) > f W(Dy) dx. 
e ^° Jn 



(2.5) 



Finally, convergence of the external force term is a consequence of Lemma 5.3 in 



BDMG99 . This result implies that if y £ — > y in L (f2) and ||.Dy e || 2 is uniformly bounded, 



then we have that 

eE £ (y £ )^ [ f( y -L(x+l))dx, (2.6) 

Jn 

and combining (2.5) and (2.6) proves Proposition [2] 

2.2. The limsup inequality. Now we have obtained the liminf inequality, we need to 
prove the limsup inequality to complete the proof of Theorem [TJ 

Proposition 4. For any y G A(L), there exists a sequence of y £ G A £ (L) such that 
y £ — )■ y in L 2 and 

limsup .F% £ ) < Foiy). 
The construction of the sequence y £ requires a diagonal argument which is similar 



in flavour to that employed in the proof of Theorem 4.5 in Bra02 . This argument 
proceeds in two steps. In the first step, the convexity of W is exploited to show that 
a naive approximation of y G A(L) by T £ y works for the 'pure' part of the energy. By 
linearising deformations near the defect, and therefore controlling the behaviour of the 
energy there, we can take a diagonal sequence to arrive at the correct inequality. Note 
that the inequality is trivial if J"o(y) = +00, so we only need consider y G A(L) such that 
T (y) < +00. 
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Fix y G A(L), and define the integrand 



N-l 

W e (x) := 

i=-N 



J2 (lMD2yi-i) + MDm) + l<h( D 2Vi)) ■ Xi(x), 

i=-N 

where Xi( x ) * s the indicator function for the interval (xi, Xj+i). Note that by construction, 

fw £ (x)dx = eE;(y). 
Jci 

We will apply Fatou's lemma to the functions W £ . 

Almost everywhere convergence ofW e . For any x G Q, define the sequence i £ := \_Nx\ G 
Z. Then 



as e — > 0, and applying Lebesgue's Differentiation Theorem (see for example Corollary 2 
in Section 1.7 of EG92| ) gives that for almost every xefl, 



DjVie = -f Dy(t) dt ^ Dy(x) 

J x 6 



as e — > 0. An immediate consequence of this and the continuity of the potentials (pi is 
that 

W £ (x) -> W(Dy(x)) 

for almost every x G Q. 

Pointwise upper bound on W £ . Fix a point i6(l and the sequence i £ as above. Dropping 
the subscript, we estimate 

W £ (x) < l<h{D 2 yi-i) + ^i(Am-i) + MDiVi) + \<h.{p 2Vi ) + IfciDzVi) + C> 
= ±W(£> 2 2/i-i) + ^(Dij/0 + §W(£> 2 Jfe) + C, 



where — C G R is a lower bound for 0i. Next, Assumption (4) in Section 1.5 implies that 
for some CgI, 

Hence, letting A := max{|, j^}, 

W £ (x) < A (w(D 2 y^ 1 ) + W(D iy A + W(D 2 yS) + C, 

i W(Dy) dx + 4- W(Dy) dx + i W(Dy) dx 1 + C 

J Xi-l J Xi J Xi J 

by using Jensen's inequality. Since W is bounded below, we can extend the domain of 
integration in each interval, possibly changing the constant C, to reach the upper bound 

rx+2e 

W £ (x)<g £ (x):=3Af W(Dy)dt + C. 

J x-2e 
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Convergence of g e . As J-a{y) is bounded, W(Dy) is integrable and so Lebesgue's Differ- 
entiation theorem implies that 

g £ (x)^g(x):=3AW(Dy(x)) + C 

almost everywhere as e — > 0. Furthermore, g e is in fact the convolution 

9e = 9*Pe, 

where p £ is an approximation to a Dirac mass given by 

Pe{x) ■= —pi —J, where p(x) = Xn(x). 

The functions g £ are in L 1 (fi) because g G L x (f2) and p £ G L°°(f2), and by standard 
arguments 



/ g £ (x)dx -> / #(x)dx 



as e — >• 0. We can now apply Fatou's Lemma to g £ — W £ , which is positive, measurable, 
and converges almost everywhere as in to get 

gdx — lim sup / W £ dx = lim inf g £ — W £ dx > / g — W(Dy) dx. 
Jn e^o Jn £ ^° Jn Jn 

A rearrangement of this inequality allows us to conclude that 

lim sup eEl (y e ) = lim sup W £ dx< / W{Dy) dx. (2.7) 

Controlling the energy near the defect. For any y G A(L) and rj > 0, let y v be a linear- 
isation close to of y G A(L) given by 

■n/ \ I y{—if) + v ™~ y }~' n > (x + rf) when \x\ < 77, 
y'yx) :— < zr > 

\y{x) otherwise. 

Let F' n := Dy v (0). For e sufficiently small, the defect energy for T £ y v is: 

eE d (T £ y^ = 2e(MF v ) - MF V ) + MF V ) - MF V )) < C(rj)e (2.8) 

with rj fixed. To control sEf, we can once again employ the estimate that was proven in 
(2.6), since the argument used was for a more general sequence than that chosen here. 
Therefore, combining (2.6), (2.7) and (2.8), we deduce that 

limsup.F £ (TV)<-W). 

Since y v — > y in L 2 as 77 — > 0, we would like to show that J-o{y v ) —^ J~o{y) as 77 — \ in 
order to use a diagonalisation argument. This follows from the observation that 

rv 
2r)-C< 



n rv 

W{Dy n ) dx< W(Dy) dx, 
n J -v 



since W is bounded below and convex. Both sides tend to as 77 — > 0, so we can deduce 
that 



W) = f WW) + f{f - L(x + §)) dx, 



12 \\i> y 112 ' 



< W(Dy)dx+ W(Dy^)+ \ f udx +\\f\\ 2 \\y* - y\ 

Jfl\(—ri,ri) J —r\ 

-> My) 

as i] — > 0, recalling that u = y — L(x + |). 
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Conclusion of the argument. Finally, by taking a diagonal sequence from the collection 
of T £ y v , there exists T e y qe — > y in L 2 , along which 

limsup^ £ (rV E )<^o(y), 
proving Proposition [4], and therefore concluding the proof of Theorem [TJ 

Remark 3. The defect does not introduce a perturbation to the r-limit at this order 



- see Theorem 3.2 in BC07 for the r-limit of this problem without a defect. This is to 
be expected, since the 'defect set' is null in the limit as e — > 0, and it therefore becomes 
reasonable to ask whether there is a higher order change in the energy, which is the subject 
of the subsequent analysis. □ 



3. Properties of F, 



o 

The functional Tq is of a well-studied form, and the analysis of the minimum problem 
is classical. The following theorem collects relevant results regarding the functional and 
its minimisers which we will invoke in the following sections. 

Proposition 5 (Properties of th -order limit). The problem 

argminj r (?/) 
veA(L) 

has a unique solution y, which has the following properties: 

(1) y satisfies the Euler-Lagrange Equations for this problem, 
(2)yeC 2 (Q). 



Proof. The existence part of this proof is completely classical, and can be found in Dac08 
for example. If we suppose for the moment that minimisers are in W 1,00 (f2) and satisfy 
the condition 

Dy(x) > 5 > 

for almost every x 6 f2, it is also easy to show that they satisfy 

f [W'(Dy) -a]Dvdx = Vt) e Wj ,00 (fi), (3.1) 

Jn 

where a(x) := J f(t) dt. Since W is /-convex, W : 1R + — > K is strictly increasing and is 
a C 1 diffeomorphism. If (IV) -1 is the inverse of W, it is possible to 'explicitly' define a 
solution of the Euler-Lagrange equations 



where £ is the solution of the following implicit equation: 

[(W')- 1 (cj(t) + Z)dt = L, 
Jn 

We can show that this equation has a solution by regarding the left hand side as a function 
of S, showing it is C 1 , has a strictly positive derivative, and tends to as £ — > — oo, so 
attains all possible values L > only once. It is now simple to verify that y satisfies the 
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Euler-Lagrange equation pointwise and is C 2 , so all that remains to do is show that this 
is in fact the minimiser. Suppose that y G A(L) minimises Fq an d is not equal to y; then 



0> / W(Dy)-W(Dy)-f(y-y)dx, 

Jn 

= [ W{Dy) - W(Dy) + a(Dy - Dy) dx, 
Jn 

> [ [W'(Dy) - a -Y,]-(Dy- Dy) dx + \l \\Dy 
Jn 



Dyf 



2 ■ 



where we have integrated by parts on the second line, and used the fact that W is /-convex 
on the last line. Since y has been constructed to solve the Euler-Lagrange equation 
pointwise, the integrand vanishes, and hence y = y. This argument clearly also implies 
uniqueness of solutions. □ 



4. The 1 st -order T-limit 



The approach taken in Section 2.2 gives a strong indication of the scaling of the next 



term in an asymptotic expansion of the energy: (2.8) suggests that the extra energy from 



the defect is only coming from a set near the defect that is of size 0(e). Section |3| shows 
that we have a very clear understanding of the properties of y, and thus we can reasonably 
hope to derive a good characterisation of the next order limit, as in |SSZ11 BC07| . 
For this purpose, we define some additional notation. Recalling from Section [3] that 

y = argminJ r (?/), 

yeA(L) 

the functional from which we obtain the first-order limit is 

To make the notation used in this section more concise, we let F := Dy(0) as in Section 



L4J and define potentials 

$!(£) := <MF +t)- MFo) ~ 4>'i(Fo) • *, 
*i(t) := MFo + t)- 0i(F o ) - ^(Fo) • t, 
$ 2 (t) ;= 2 (F O +t)- MFo) - MFo) ■ t, 
*a(*) := MFo + t)- MFo) ~ MFo) • t 
We will show that the first-order r-limit can be written in terms of the infinite cell problem 

inf E^(r), 



where E n 



(ZW 



re£ 2 {1) 

U {+c»} is defined to be 



E n 



E ^ 



$2 



"Ti + l+Ti 



)+E d ( 



and we have set 
E d {r) := iJ 2 (F + 



r-2+r- 



l ) ~ M F o + r -^Y^) + Vi {Fo + r_i) - 0i (F + r_x) 



+ fa (F + r ) - 0i (F + r ) + fa (F + 



rp+ri ' 
2 



; (F + ^±a). 
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The second main result of this paper is the following theorem. 

Theorem 6 (l st -order r-limit). With respect to convergence in L 2 , we have that 

( inf -Eoo(r) y = y, 
r-\JmFZ(y)=r 1 (y):={ret'W 

£ ^° l+oo y^y. 



In contrast to the results of SSZ11, BC07 , we emphasise that we have an explicit 



representation of the l st -order limit in terms of a minimisation problem in an infinite 
cell. Once more, the proof of this result divides into two parts, the liminf and limsup 
inequalities, which we prove in the next two sections. 

4.1. The liminf inequality. The liminf inequality is the following statement. 

Proposition 7. If y e —¥ y in L 2 , then 

liminf 7^ £ ) > Ti(y). 

As in the proof of Proposition [2j we use a coercivity result which says uniform bounded- 
ness of J-'Ky 6 ) implies a form of compactness. In this proof, there are two such results, 
which are employed at crucial steps in the main argument. The first of these results, 
Lemma [8j states that if J-'fiy 6 ) is uniformly bounded then the weak convergence of y £ in 



LI 1 proven in Lemma [3] improves to strong convergence in H 1 . The second, Lemma 10 
describes coercivity in a topology which we use to describe perturbations to the minimiser 
of J-o close to the defect. Once these results have been obtained, the main argument will 
follow by applying Fatou's Lemma to a suitable reinterpretation of T\ (y £ )- 

The first key step before proving the coercivity results is to rewrite J-|(y £ ) by using 
integration by parts on the external force terms. For J~o(y), 

[ f(y-L(x + I)) dx = - [ a(Dy - L) dx, 
Jn Jn 



using the boundary conditions, where a(x) := f_ x / 2 f(t) dt as in Section 1.4 Analogously, 
recursively define 

0~-N := —2 e f~Ni 

This leads to the representation 



JV-l N-l _ _ N-l 

CT7, i — 0"f 



e 

i=~N i=-N i=-N 



N-l 



i=-N i=-N 

We define the step function a £ : Q — > 

° £ {x) ■= Yl a !+iX(xi,x i+1 )(x), 

i=-N 

so that if y e A £ {L), 

eEI{y) = - I a £ Dudx. 
Jn 
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Using these definitions, we perform careful estimates of T\ (y £ ) by splitting the domain of 
integration over the intervals (xj, 2^+2) • For y £ G A £ (L), define 



J Xi 



s £ := HDiVf) + iMDiy!) + lMDiyf +1 ) - \o £ +1 D lU £ - \a £ l+2 D lU £ +1 

W(Dy) -aDu + T, (Dy £ - Dy) dx. (4.1) 

if i G {— N, . . . ,N — 1}, and set s £ := otherwise. Then it is easy to check that 

00 

■W) = Y, s£ i + M D iy £ - F °)- ( 4 - 2 ) 

i=— 00 

We are now in a position to prove the coercivity results. 

Lemma 8 (Strong H 1 coercivity). If y £ — > y in L 2 , and J^fiy 6 ) is uniformly bounded 
for all e > 0, then y e — > y in H 1 . 

Proof. Since Ffiy 6 ) is uniformly bounded, we know that for some CgI, 

\jr%tf)-My)\<Ce, 

which immediately implies that ^(y 6 ) — > J-o(y) as e — > 0. Consequently, Lemma [3] 
applies and so ||-Dy e || 2 is uniformly bounded. Let i G {— N, . . . , N — 1}. We estimate sf 
below: 



/ ' 2 (w{Dy £ ) - W{Dy) - a £ Du £ + oDu - £ (ity 6 - £>y)) dx, 

/ 2 (V'(£>£) (Dy 6 -Dy)+l \Dy £ - Dy\ 2 - a £ Du £ +aDu-E (Dy £ - Dy)) dx, 

J Xi 

using the concavity of 2 on the first line, and the /-convexity of W on the second. Next, 



s £ > 

' x 

'■X 

> 

' X 



the Euler-Lagrange equation (|3.1|) implies that 

s £ > 



fXi+2 . . 

+ f 0"(Dy £ -Dy) +1 \Dy £ - Dy\ 2 - a £ Du £ + cxDwj dx, 

J Xi 

>-f l2 ((a-a £ )Du £ + l \Dy £ - Dyf\ dx. (4.3) 

J Xi 

The latter term in the above integral is of the form we are looking for, so it now remains 
to show that the other term vanishes in the limit. Once this is done, we then show that 
the defect energy is also suitably bounded below. 

Pointwise estimate on a — o~ £ . Noting that Du £ is constant on the intervals (xi, JCj+i) and 
using the definitions of a and o~ £ , we rewrite 

/ ! (a-a £ )dx = j ' ( / /(()dt-ej Wi + fj+i) ~ ¥h ) dx, 

t/x, ix, \ J —1/2 j=-N ' 

l i2 (f(t)-{T £ f)(t))dt+j- (J f{t)dt-\efAdx. 

Since we know that / G C 2 , standard results about interpolation error (see for example 
SM03| ) imply that 



f[ 2 {f{t)-{T £ f){t))dt 



<^nriL- (4-4) 
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For the other term, we Taylor expand f(t) at Xi, then evaluate integrals to show that 

"X i+1 rx 



rxi+i rx 

f / f(t)dt-lefidx 

J Xi J Xi 



< -e 2 \\f'\\ 

— 6 II -/ lloo 



(4.5) 



Combining (4.4) and (4.5), we have 



J Xi 



(a — a £ )Du £ dx 



<(^ 2 \\f"\L + ^ 2 \\f\ 



\D lU f\ + \D lU f +1 \ 



<Ce 2 (\D 1 u £ \ 2 + \D 1 u 
< Ce 3/2 \\Du £ 



e |2 
j+ll 



1/2 



l 2 , (4-6) 

where on the second line we used Jensen's inequality, and on the third line we used e 1 ^ 2 
and added further postive terms inside the brackets to get the estimate. This can be used 
in (4.3) to give 



s £ > -Ce 3/2 \\Du% + 1 1 2 \Dy £ - Dy\ 2 dx. 

J Xi 



(4.7) 



Lower bound on defect energy. Using (2.2), the fact that W(Dy) — (a + T,)Du is finite 



and estimate (4.6), 

o 



J2 at + E d (Dy £ ) > J2 \ l {\ D ^\ 2 + l^iyf+il 2 ) ■■"<" 



i=~2 



i=-2 





Xi+2 



> 



a £ Du £ + ZDu £ )dx, 

J2 -f 2 ( i \ D v e \ 2 - aDu£ - (° £ - °) Du " - s DuE ) dx + °> 

i=-2 ^ 

J2 f ! 2 (l \Dy £ \ 2 ~ Ik + EIL \Du'\ + C)dx- Ce 3 ' 2 \\Dy% . 



i=-2 Jx 

Since Dy is bounded above and below, by adjusting constants suitably we have that 
o o 



f*Wn I 9 

Y, 4 + E d (Dy £ ) >J2f l \°y £ -Dy\ 2 dx + C- Ce 3 ' 2 \\Dy £ 

i=-2 i=-2-> x i 



|| 2 . (4.8) 



Conclusion of the argument. By summing over i in (4.7), combining with (4.8), and using 
the fact that ||Z}?/ £ || 2 is uniformly bounded, we have shown that 

X- 1 fx i+1 

H{y £ ) > -Ce 1 ' 2 + i Y,{ W- D v\ 2 dx + a 



-N° Xi 



Multiplying this inequality by e and using the assumption that J~[(y £ ) is uniformly 
bounded, we have 



— II 2 



Ce>l\\Dy £ -Dy\\ 



2 ■ 



which proves the result. 



(4.9) 

□ 



To prove the second coercivity result, we define the sequence of operators V £ : A £ (L) — > 
? 2 (Z)by 



(P<v)i - 



D lVi -D x y { ie{-N,...,N-l} 







otherwise. 
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Clearly V e y is well-defined since this sequence is non-zero only on a finite set. 

Lemma 9 (Weak £ 2 (Z) coercivity). If J^Ky 6 ) is uniformly bounded, then there exists 
a subsequence ofV £ y £ which converges weakly in £ 2 (Z). 



C>lJ2i \Dy e - Dy\ 2 dx > I J2 f Dy £ -Dydx 



Proof. By dividing (|4.9|) by e and using Jensen's inequality, we have 

-•N ox i i=-N J Xi 

We have shown that the sequence V e y e is uniformly bounded in £ 2 (Z), so in particular, it 
must have a weakly convergent subsequence. □ 

To conclude the argument which will prove the liminf inequality, we will use the fol- 
lowing characterisation of weak convergence in £ 2 (Z) which follows easily from the Riesz 
Representation Theorem. 

Lemma 10. A sequence (r £ ) C £ 2 (Z) converges weakly to r G £ 2 (Z) as e — >■ if and 
only if the following two conditions hold: 

(1) 1 1 r £ 1 1 £2 is uniformly bounded, 

(2) r e —y r pointwise (almost everywhere in the counting measure) as e — y 0. 



As indicated at the beginning of this section, we apply Fatou's Lemma to the sum (4.2). 
Suppose that J-'Ky 6 ) is uniformly bounded and y e — >■ y. Take a subsequence y £k such that 

lim JFl h {y eh ) = liminf JR[{y E ), 

k— >oo e— s>0 

and then using Lemma [9j a further subsequence (which we do not relabel) such that 
V £k y £k weakly converges to r in £ 2 (Z). Since y G C 2 (f2), we have that 



D 1 y i = f " Dydx^F :=Dy(0) 

J Xi 



as e — y 0. Fixing an index ieZ, Lemma 10 implies that 

(V £k y £ *) i = D 1 y £k -D 1 y l , 

^Fo + n-F 

as k — y oo, so that we may view r as a perturbation to the deformation gradient in an 
'infinitesimal' neighbourhood of the defect. The 'pointwise' estimate (4.7) implies that 
forie{-JV,...,iV-l} 

s\ + Ce 3/2 > 0, 

so that 

A/fe-l oo 

liminf V s £k + Ce 3 k /2 > V liminf (s £k + Cel /2 ). (4.10) 

i=—N k i=—oo 

Since the potentials 0j are continuous and erf — )■ cr(0) as £ — >■ with % fixed, we have that 
liminf (s £k + ef) = 2 (F O + ^p±i) + §0^0 + r<) + §0i(F o + r i+1 ) - W(F„) 

- cr(0)(F + ^±i) + <r(0)F - E ^f±i. 
Recalling that y satisfies the Euler-Lagrange equations pointwise, we have that 

<t(0) = W(F ) - E, 
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so that 

liminf « fc + sf) = $ 2 (r, + r l+1 ) + i^fa) + |$i(r m ). (4.11) 

Note that 

-3/2 _ ,- ~ 1/2 



lim V CeJ /2 = lim C £ * /2 = 0, (4.12) 

fc— >oo ' ^ fc— >oo 

i=-JV fc 

so then combining (4.10[), (|4.11[) and (|4.12), we have that 



liminf J2 s ? > J2 M 1 ^) + iMn) + l$i(r i+1 ). 

i=—oo i=—oo 

Finally, by possibly taking further subsequences, we can assume that (V ek y Sk ) . converges 
uniformly for i G {— 2,— 1,0,1}, and then we have 

liminf Fl{y £ ) = liminf ( V s^ + E d {Dy Ek ) j , 

^ i=— oo 7 

> ^oo(r), 

> inf -Eoo(r), 

proving Proposition [7} 

Remark 4- By definition, we have 

2V-1 

E (^ £ ), = ° 

i=-JV 

for any y £ G *4 e (L). If P e y e — *■ r in ^(Z), we could conclude that 



I> = 0; 



2= — OO 

02/ 



however, since we have convergence only in £ (Z), this is not true in general. It is therefore 
clear that the set of compactly supported mean zero sequences is dense in the weak 
topology on £ 2 (Z). □ 

4.2. The limsup inequality. The limsup inequality is the following statement. 

Proposition 11. For every y G A(L), there exists a sequence y £ — >■ y in L 2 such that 

lim sup J={{y e ) < F x (y). 

This statement is trivial in the case where y ^ y, so we only need to construct the 
sequence for y — y. In order to construct the limsup sequence, we will show that there 
exists a minimiser of Eoo, and then combine a suitable truncation of this minimiser with 
T e y to get the result. 

Proposition 12 (Properties of l st -order limit). Let E^ : £ 2 (Z) — > R be defined 
as in Section y. Then there exist minimisers of E^ which satisfy an infinite system of 



nonlinear algebraic Euler-Lagrange equations, given in (4.13). 
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Proof. Since i?oo( r ) < +00 for the constant sequence r = 0, the infimum is less than +00. 
We show that existence follows from the direct method of the Calculus of Variations 
applied to E^. It is easy to check that $2 is concave because 02 is concave, hence 

00 

^oo(r) > Yl $ i( r *) + $ 2(r.) + E d (r). 

i=— 00 

Next, 

&i(t) + $ 2(t) = W(F + t)- W(F ) - W'(F )t > \lt\ 



using the /-convexity of W . For i G {—2, —1, 0, 1}, we can estimate below as in (2.2), so 
that for some constant C G M. 

1 1 

Y, (§*i(r«) + i*i(r m ) + $ 2 (^)) + E d (r) > ^ \l \r t \ 2 + C. 

i=-1 i=-2 

Hence we have that 

Eoo\ r ) > 2 H r ll^ 2 (z) "I" ^- 

Next we need to show that E^ is sequentially weakly lower semicontinuous. This follows 
by using Fatou's lemma as above with the pointwise lower bounds just proven. A standard 
application of the direct method now yields existence. 

To obtain the Euler-Lagrange equations, suppose r is a minimiser of E^. Let e l G £ 2 (Z) 
be the sequence which has 

3 1 otherwise. 
Let i <£ {—2, —1,0, 1}; for small enough t > 0, E^r + te 1 ) < +00, and 
Q< Epo(r + te«)-.E 00 (r) 



t 

= [ l& 2 ( ri ~ 1+ 2 n+at ) + $ifa + st) + l& 2 ( ri+s T i+1 ) d«- 
Jo 

Applying the Dominated convergence theorem and repeating the argument for t < now 
implies that 

1^(2=^) + $i(r0 + f<& 2 (^) = 0. (4.13a) 

By the same argument, we also have that 

i& 2 (!=*£=*) + $;(r_ 2 ) + \%{ T -^) = 0, (4.13b) 

|^ 2 (!^±-i) + ^ r _ x ) + |^(^fi) = 0, (4.13c) 

i$ 2 (^o) + ^(ro) + |*' 2 (^) = 0, (4.13d) 

iv]/' 2 (ro±n) + $; (ri ) + |^(n±a) = , (4.13e) 

completing the proof. □ 



In order to complete the proof of the limsup inequality, we will require a better under- 
standing of minimisers of E^, and so we prove the following sequence of results, which 
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amount to regularity results for solutions of the Euler-Lagrange equations (4.13). From 
now on, we fix r G £ 2 (Z) as being one particular minimiser of E^. 



Lemma 13. Suppose that r G £ (Z) solves (4.13a) for all i > 2. Then 



r\ > =>- ri = max{r,}; 

jeN 

n < =>• n = min{r.,j. 

Proof. We prove only the first conclusion, the proof of the second being similar. Suppose 
for a contradiction that tm is an interior maximum, i.e. that 

Tm = maxjr,,} > 0. 
Then because $2 is concave, we have that $ 2 is monotone decreasing, and hence 

= I^TM^IM.) + $' l(rM ) + i$' 2 (^±^), 

> $ 2 (r M ) + $i(nf), 

= ^'(F + r M )-^'(F ), 

iy"(F + s) ds, 


>/r M , (4.14) 

which implies that either rj = for alH > 1 or a contradiction, concluding the proof. □ 

Corollary 14. If ' (r,)^ 1 G £ 2 (N) solves the Euler-Lagrange equations with r x /ixed, then 

ri > ^> rj > fj+i > /or all i G N; 
ri < =^ ri< r i+ i < for all i G N. 



Proof. Estimate (4.14) states that if 

r i = m&x{r i _ 1 ,r i ,r i+1 }, then n < 0. 
Similarly, it is possible to show that if 

ri = min{rj_i, r h r i+i }, then r t > 0. 

Suppose that r has a strict local maximum ru < with M > 1. Then r^+i < rjif < 0. 
Since local minima can only occur when r^ > 0, rM+i cannot be a local minimum, and so 
fM+2 < fM+i < 0. Proceeding by induction, Tj < r^ < for all i > M, which contradicts 
the fact that r G £ 2 (N). A similar argument prevents the existence of strict local minima. 

Next suppose that tu = is a local maximum for M > 1. If Tm+i < 0, then the 
previous argument applies. If r^+i = 0, then it too must be a local maximum or minimum, 
depending on the sign of r M+2 . If r M+2 7^ 0, then we can apply the previous arguments 
again to arrive at a contradiction, so by induction we have that r, = for all % > M. 

We have therefore shown that there can be no internal maxima, unless they are degen- 
erate in the sense that r is identically after the maximum, and by a similar argument, we 
can show that there can be no internal minima except if they are degenerate in the same 



sense. We can now conclude that any solution of (4.13a) must be increasing if n > 0, or 



decreasing if r\ < 0, which concludes the proof. □ 
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Finally, we prove that minimisers have exponentially small 'tails'. 

Proposition 15 (Exponential decay). Let C > be the constant such that for all 
te (-ri,ri), 



> $' 2 '(i) > -C. 
i+c> 



Then if X :— -A;, we have that 



W < X l ~ In 



Proof. We will only prove the result for n > 0, since the other case is similar. The 
Euler-Lagrange equations may be rewritten using the Fundamental Theorem of Calculus 

as 

f l W"(F Q +t)dt+ [ ' $2(t) dt+ f $2(t) dt = 0. 
Corollary |14| now gives that r^-i > r, > r^+i, so we have 

= f % W"(F + t)dt+ f l $%(t) dt- J $%(t) dt, 

Jo J n Jr i+ i 



'r i+ i 

> I i dt- i ' l C dt. 



using the assumed bound on the second derivative of $ 2 , and the Z-convexity of W. It 
immediately follows that 

C 



l + C 

which is true for any i > 2 and the decay estimate is obtained by using induction on this 
inequality. □ 



Remark 5. It should immediately be noted that since we know that r, converges to 
zero exponentially as i — > ±oo, it must be the case that r G £ 1 (Z). This will be crucial in 
what follows. □ 

We can now apply this characterisation of minimisers of E^ to complete the proof of 



Proposition 11 Let the sequence of functions y e ' s be given by 

y £ '\x):= [ X Dy*> s (t)dt, 

J -1/2 



-1/2 

where we have set 



DxHi otherwise, 

K-l 

f s -.= 6j2 n, 

i=-K 
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We will prove estimates for y £,s , and then set 6 as a function of e in order to obtain the 
recovery sequence. Since r e -^"(Z), we have that 



K-\ 



\f s \ < 



8 2_, \ r i\ — $ \\ r \ 



^(Z) 



(4.15) 



i=-K 



By construction, y £)S £ -A £ (L) as long as 5 > e, i.e. K < N. 



Ensuring y £,s is well-defined. The first step we take is to check that W{D\yl' ) is well- 
defined. Since E^r) < +oo, we have that for all i G Z, 

W(F + n) < +oo =* F + r* > F - ||r||^ (z) > 0, 



as £ 2 (Z) C £°°(Z). Next, by using the Mean Value Theorem and (4.15), we estimate for 
i e {-K, . . . , K - 1} that 



Dxyf -F -n 



Zi+i 



Dydx-r 6 -Dy(0) 



< 



\Dy(Z)-f s -Dy(0)l 
D 2 y(x) dx 



+ S\\r\ 



e 1 (z) 



< \x 



i+l\ 



D'y\L + 5\\r\\ el 



(Z) 



< ^ Felloe + 5 H^ 



^(Z) • 



On the second line, £ is some point in (xi,Xi+i), and on the final line we have used 
\xi + \\ < Ke = £j. For fixed i, it is now clear that 

\W{D lV f) - W(F + n )\ < Csup {\W'(t)\ :te(F + r h D l3 /|)}(| + <J), 

therefore W(Diyf) is finite when e<(J<l. We can additionally estimate 



iDij/^-Dj/Or) 



/•asi+i ,• v 

-/- ( Dy dx — r s — Dy{x) j da: 

•/Xj 



<e||D 2 v|| + 5IH 



P(Z) ■ 



(4.16) 



^s,S 



Pointwise upper bounds on s^ . Next, we define s t ' in a similar fashion to (4.1), but 
adding and subtracting an extra (^xiD^yl' ) term, we have 



„e,<5 / xxr/ n „.£:<5 



,; = \W(D 2 yf°) 



Xi+2 



W{Dy) dx 
+ -[ 2 (oDu - o- £ Du £ ' 5 - £ (Dy £ > s - Dy)) dx, 

JXi 



=:Ti+t 2 +t 3 . 
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T\ can be treated using Jensen's inequality and the convexity of W: 



T 1 = W(D 2 y £ '°)-f W(Dy)dx 

JXi 

< \W(D iy f) + \W(D x y^ x ) - f ,J W(Dy) dx, 

JXi 

< I 2 W{Dy £ > & ) - W{Dy) dx, 

J Xi 

< I X W\Dy £ ' S ) (Dy £ ' s - Dy) dx. 

J Xi 



(4.17) 



The final integral term, T 3 , is bounded by 



To 



%i+2 



aDu - a £ Du £ ' 6 - £ (Du £ ' 5 - Du) dx 

i 

I "V + s ) (Du - Du £ ' 5 ) + (a- a £ )Du £ ' S dx 

J Xi 

I 2 W' (Dy) (Dy - Dy £ ' 5 ) dx + Ce 3/2 \ \ Du £ ' s \ 

J Xi 



(4.18) 



using the Euler-Lagrange equations and the estimate proved in (4.6). The remaining part 
of the expression, T 2 , can then be estimated as follows: 

t 2 = iMDiyf) + \MDwtt) - M D *yf) 

v;( s -¥)dsdt, 



J Diy i ' JD lVi ' 

< C IDiy^ - D lVi ' | , 

where we have used the fact that </>'/ is bounded on the domain of integration for sufficiently 
small e and 5, and (f)\ G C 2 . If i G {— K, . . . , K — 1}, then applying the Mean Value 
Theorem with £ G (xj,Xj + i) and ( G (x i+ i,x i+2 ) gives 



2 



't+l| 



^i(Al/f) + iMDwtti) ~ M D 2Vi) < C\Dy(0 ~ Dy(0 + r i+1 - 

<C r Q|D 2 y|| oo e+|ri| + |r i+1 
<cf ||£> 2 y|| 2 e 2 + |r. 

— V II Moo ' 

In the case where % ^ {— K, . . . , K — 1}, we obtain 

Ifcpitf') + ^ipiygi) - H D *yf) < c z 2 \\ D2 y\L ■ 

Adding together J4.17| ), ( |4.18[ ) and ( |4.19[ ), for i G {-K - 1, . . . , K - 1} we have 
/ 2 (V(Zty £ ' 5 ) - WW)) (W'* - ^y) dx + Ce 3/2 \\Du £ > s \\ 2 

J Xi 

+ Ce 2 + C(\r t \ 2 + \r i+l \ 2 ), 
C-F** \Dy £ > 5 - Dy\ 2 dx + c(e 3/2 + e 2 + \r^ + 

<c[{s + 5) 2 + e 3 / 2 + e 2 + \r i \ 2 + \r i+1 
<C7(5 2 + e 3 / 2 + e 2 + |r i | 2 +|r. 



(4.19) 



sf < 



< 



R+il 



»+i 
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using W G C 1 , estimate (4.16) and Jensen's inequality. For the other indices, we obtain 

sf < I 2 (w'(Dy £ ' s ) - W'{Dy)\ (Dy £ ' S - Dy) dx + Cs 3/2 \\Du £ ' s \\ 2 , 

<C(e 2 + s 3 / 2 ). 
Conclusion of the argument. The two pointwise estimates just obtained imply 

oo JV-1 K 

J2sf< C J2^ /2 + Ci: (\r,\ 2 + 5 2 ), 

i=—oo i=—N i=—K 

<c{e l / 2 + \\r\\l m + 5). 

By choosing K = \_yN\ , we have that S < Ce 1 ^ 2 , and so an application of Fatou's Lemma 
with the pointwise upper bound we have just proven implies that 

oo oo 

limsup y sf + Ed(Dy £ ) < > limsups^' + \imEd(Dy £ ), 

l=— OO l = — oo 

= L(r). 



This now proves Proposition 11, and concludes the proof of Theorem pi 



Remark 6. This result shows that the perturbation to the minimiser from the con- 
tinuum model is confined to an exponenti ally thin boundary layer. Note that the lin- 



earisation of the functional E^ in Section 1.4 yielded a similar solution structure; this 



exponential decay suggests that any interaction between defects of the type described 
here is likely to 'decouple' if one were to study a situation in which there were multiple 
defects of a fixed and finite number which are well-separated in the limit e — > 0. □ 

Conclusion 

We have presented an analysis of a model for a point defect in a ID chain of atoms 
interacting under assumptions which attempt to replicate a Lennard- Jones type interac- 
tions in an elastic regime. We have derived the th -order r-limit, which is identical to the 
limit when there is no defect. 

We then proved that the l st -order r-limit exists and have given an explicit character- 
isation of this limit in terms of an infinite cell problem, and shown that the perturbation 
introduced by the defect is confined to an exponentially thin boundary layer. 
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